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A TWO-DIMENSIONAL OBLIQUE EXTENSION OF BESSEL 

PROCESSES 


Dominique LepingkQ 


Abstract. We consider a Brownian motion forced to stay in the quadrant by an electro¬ 
static oblique repulsion from the sides. We tackle the question of hitting the corner or 
an edge, and find product-form stationary measures under a certain condition, which is 
reminiscent of the skew-symmetry condition for a reflected Brownian motion. 


1. Introduction 


In the present paper we study existence and properties of a new process with values in the 
nonnegative quadrant S = M+ x M + where R + := [0, oo). It may be seen as a two-dimensional 
extension of a usual Bessel process. It is a two-dimensional Brownian motion forced to stay in 
the quadrant by electrostatic repulsive forces, in the same way as in the one-dimensional case 
where a Brownian motion which is prevented from becoming negative by an electrostatic 
drift becomes a Bessel process. Note here and now that the corner 0 = (0,0) will play a 
crucial role and in some cases it will be necessary to restrict the state space to the punctured 
nonnegative quadrant S° = S \ { 0 }. 

Let (fi, T, (Pt)t> 0 i P) be a complete probability space endowed with a filtration (Ft)t >o 
satisfying the usual conditions. Let (B tl Ct) be an adapted driftless Brownian motion in the 
plane starting from 0, with covariance matrix 

1 P 
P * 1 

and p € [—1, +1]. 


Definition 1. Let a, /?, 7 , h be four real constants with a > 0,(5 > 0. We say that an (Pt)- 
adapted continuous process (X,Y) with values in S is an Oblique Two-dimensional Bessel 
Process (02BP) if for any t> 0 

m X t - X 0 + B t + afi ^ + & fa y- e > 0 

Y > = r„ + c t + 7 / 0 'f+ > 0 

where A"o and Yq are nonnegative J~o-measurable random variables, and 

Joj 1{A S =0 }ds = 0 ./o l { y s=0 }ds = 0 

Jo -*-{A s >0}A7 < 00 Jo 1 {U>0}^ < 00 ■ 


This stochastic differential system is very singular at the edges of the quadrant and the 
question of existence and uniqueness of a solution is not simple. The particular case when 

b = 7 = 0 

(2) U t = Xq + B t + a fg jf- s > 0 

V t = Y 0 + C t + 6f*^ > 0. 
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is already known: the processes U and V are Bessel processes. Actually, U is a Bessel process 
of dimension 2a + 1, and the point 0 is instantaneously reflecting for U if a < \ and polar if 
a > i. If p = 0, U 2 + V 2 is the square of a Bessel process of dimension 2a + 25 + 2 [32], and 
so the corner 0 is polar for (U, V) in this case. Comparison between X and U, Y and V will 
play a key role in the construction of the solution (X, Y ) and the study of its behavior close to 
the edges of the quadrant. The process (U, V) is an example of Brownian motion perturbed 
by a drift deriving from a convex potential. More generally, stochastic differential systems 
including such a singular drift have been studied in [2, 26, MIEZ], where strong existence and 
uniqueness were obtained. They are examples of so-called multivalued stochastic differential 
equations, also called stochastic variational inequalities in convex analysis. 

We will use these results to study the solutions to (JT]) in the oblique case where /? or 7 
do not vanish. We obtain strong existence and uniqueness for a large set of parameters and 
initial conditions, but not for all possible values. In the proofs we naturally fall into the 
crucial question of hitting the corner, that is the non-smooth part of the boundary. Using 
McKean’s argument on the asymptotic behavior of continuous local martingales obtained by 
time change from the real driftless Brownian motion, we are able to state several sufficient 
conditions to prevent the processes from hitting the corner. Our methods are not powerful 
enough to allow for necessary conditions. However we shall not restrict to processes avoiding 
the corner and, depending on the parameters, we will get existence and uniqueness (in a 
strong sense) sometimes in the whole quadrant, sometimes in the punctured quadrant, that 
is, the quadrant without a corner. 

We will also obtain some partial results about the attainibility of the edges of the quadrant. 
It is interesting to see whether the boundary behavior of one component may be modified by 
the interaction with the other component. 

In the one-dimensional case, the so-called scale functions transform the solution of a sto¬ 
chastic differential equation into local martingales and one may infer some information on 
the boundary behavior of the process. This technique was very successful in the study of 
Bessel processes. Here we still obtain functions of an 02BP which are local martingales or 
supermartingales for some values of the parameters and we derive some information on its 
asymptotic behavior. 

The laws of Bessel processe with different parameters [28] are mutually abolutely contin¬ 
uous when their paths do not reach the origin. Here we obtain two partial results about the 
absolute continuity of the laws of 02BPs for some values of the parameters. 

Finally we follow the way explored in [23], where a drifted Brownian motion is endowed with 
another drift term that is continuous and depends obliquely, via a regular potential function, 
on the position of the process relative to an orthant. Under an additional condition, which 
is called a skew-symmetry condition, an invariant density was given in an explicit product 
form. We obtain the same result for an 02BP (where the Brownian motion is now drifted) 
under an anologous condition. This time, the invariant density is the product of two gamma 
densities, which is consistent with the one-dimensional case where the invariant measure of 
a drifted Bessel process is a gamma distribution. 

1.1. Comparison with obliquely reflected Brownian motion. Studying 02BPs makes 
appear a strong connection with the properties of a semimartingale reflecting Brownian mo¬ 
tion in the quadrant. We briefly recall the definition. Let (Bt , Ct)t> 0 be the Brownian motion 
in Definition |T] 

Definition 2. An semimartingale reflecting Brownian motion (SRBM) in the quadrant is a 
continuous adapted process (. X t ,Y t ) which is a solution to the system 

/o\ Xt = Xq + Bt + L] + r\L-j: > 0 

Yt = Yq + Ct + ^2 L\ + Lf >0 
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where L 1 and L 2 are two continuous adapted nondecreasing processes with Lg = Lg = 0 such 
that for any t > 0 

[ l{x s >o }dL\ = I l{ Ys >o}dL 2 s = 0. 

Jo Jo 

Here Xq and Yq are nonnegative X^-measurable random variables, r\ and r 2 are real numbers. 

There was an extensive literary output on that topic in the eighties, with a more general 
domain : a wedge, an orthant or a convex polyhedron. We mention the works of Harrison, 
Reiman, Varadhan, Williams, Dai [3 KI3S1 SSI S] to cite a few of them. For a more complete 
bibliography we refer to [40] . 

We note a first analogy between their results and ours: we prove a necessary condition of 
existence of an 02BP is the existence of a convex combination of the directions of interaction 

(4) rx = ( 7 ) and ry = ( 5 ) 

that points into the quadrant from the corner. In the reflection setting it was proved in m 
that a necessary an sufficient condition for the existence of a SRBM is the existence of a 
convex combination of the directions of reflection 



with the same property. 

Another analogy is met in the delicate and important question of attainibility of the corner. 
The authors in [36] have found an explicit harmonic function that provides a full answer to the 
question of hitting the corner for a reflected Brownian motion in a wedge of angle £ E (0, 2 - 7 r) 
with the identity matrix as covariance matrix. We do not have any such convenient function. 
However, our second condition 

(5) 2P ~ f + a 

in Corollary [10] is reminiscent of the necessary and sufficient condition in [29] which writes 

(6) 2p < n + r 2 

in the setting of Definition [2j Restricting to £ E (0, n), we may consider a linear transforma¬ 
tion in the plane that changes the wedge in [36] into the nonnegative quadrant and the initial 
Brownian motion into a Brownian motion with covariance coefficient p = — cos£. Then the 
non-attainibility condition in [36] is transformed into the above condition ([ 6 |). 

For theoretical as well as practical reasons, a great deal of interest was taken in the question 
of recurrence of the Brownian motion with a constant drift vector and oblique reflection, and 
in the computation of the invariant measure [Mi mi Ei Ei]. Under the assumption that the 
directions of reflection satisfy a skew-symmetry condition, it was proved that the invariant 
measure has exponential product form density [TO], 11 , [39] 00]. This result has a practical 
interest because of the interpretation of the SRBM as approximation of the queue length 
processes for networks of queues in heavy traffic (7[ [10]. There is also a financial reason 
for studying SRBMs. Motivated by the so-called Atlas model of equity market presented in 
0 , some authors [f3. jj . HU El E>- 29] have recently studied competing Brownian particles 
on the line with rank dependent local characteristics. In fact the gaps between adjacent 
particles are proportional to the components of a SRBM in an orthant. There is an invariant 
probability density with an explicit exponential product form when the volatility coefficients 
are constant [24]. More generally, this is still true if their squares depend on rank linearly since 
one may infer from the statements in Section 2 of [29] that the skew-symmetry condition is 
still satisfied. Following the way in mm and [ 8 j, a more general question is the recurrence 
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or transience of an 02BP. Answering this question does not seem to be an easy task. The 
method in [12], which provides a full answer for the obliquely reflected Brownian motion 
in the quadrant, appears to break down here. Mimicking the computation in [23], we just 
calculate an invariant measure in product form under a skew-symmetry condition which is 
the equality condition in the inequality (f5j) . whereas the skew-symmetry condition in the 
reflecting case is the equality condition in ([6]). Now the terms of the product are gamma 
distributions with explicit parameters. 

Another topic of interest in [23] lies in its Remark 4.12: when a scale parameter goes to 
zero, the exponentially reflected Brownian motion should converge to the obliquely reflected 
Brownian motion. It is typically a penalty method. This kind of approximation by a sequence 
of diffusions with regular drifts living on the whole Euclidean space has been used for instance 
in [2011301J for oblique reflection in domains with smooth boundary and in [22] 33] for normal 
reflection in convex domains. With regard to our framework, it can be shown that a sequence 
of real Bessel processes with dimension decreasing to 1 pathwise decreases to a reflected real 
Brownian motion. It could be interesting to consider a sequence of oblique two-dimensional 
Bessel processes with interaction matrices 


n 


-i 


l n \ 

r 2 1 ) 


and to let n go to infinity. It should converge to the SRBM in Definition [2] This time, this 
would be an interior approximation, as it was done in m and [25] in the particular setting 
of stationary Markov processes associated with Dirichlet forms. We leave aside this point for 
further investigation. 


1.2. Organization of the paper. The rest of the paper is organized as follows. In Section 
2 we recall some trajectorial properties of usual Bessel processes. In Section 3 we state and 
prove three main lemmas of repeated use in the sequel. We also recall an existence and 
uniqueness result for a special case of multivalued stochastic differential equation that will be 
useful in our construction of an OB2P. Section 4 is devoted to the proof of sufficient conditions 
to avoid the corner of the quadrant. The main theorems of existence and uniqueness of an 
02BP are given in Section 5. In Section 6 we discuss the question of hitting the edges of the 
quadrant. In Section 7 we consider two particular cases where there exist simple functions of 
02BPs that are local martingales and obtain some information on the asymptotic behavior of 
the paths. We leave the trajectorial point of view in Section 8 to tackle questions of absolute 
continuity of the law of an 02BP with respect to the product of laws of real Bessel processes. 
The final Section 9 introduces a skew-symmetry condition that allows us to obtain existence 
of a stationary probability in form of the product of two gamma distributions. 


2. Some properties of Bessel processes 

A Bessel process of dimension d > 1, starting at r > 0, is the unique solution to the 
stochastic differential equation 

d- 1 /'* 1 

(7) R t = r + W t -t-— / — ds . 

^ JO 4t s 

where VP is a standard driftless real-valued Brownian motion starting at 0. 

We know that: 

• d > 2 : the point 0 is a polar set ([28]. Proposition V.2.7); 

• d = 2 : limsup^oo Rt = +oo, liminf^oo Rt = 0 ([28], Theorem V.2.8); 

• 1 < d < 2 : the point 0 is instantaneously reflecting ([28j. Proposition XI.1.5). 

More precisely (|T8j, p.337 and p.339), 
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• Dimension d > 2: 

P(i?t > 0,Vt > 0) = 1, 

P (Rt oo, t —>• oo) = 1, 

i? 2 ~ rf is a local martingale. 

• Dimension d = 2: 

P(i? 4 > o,vt > 0 ) = 1 , 

P(sup 4>a Rt = oo, inf t>a Rt = 0) = 1 for any a > 0, 
In R is a local martingale. 

• Dimension 1 < d < 2: 

P{Rt > 0,Vt > a) = 0 for any a > 0. 


Corollary 3. VF/ien d>2, for any a > 0, 


r°° ds 

/ -n2=°° a - S - 

J a 1L s 


Proof. From Ito’s formula 

In R t 




d —2 ft ds 
2 Ja Rj' 


From ([28]) Proposition IV.1.26) we know that on the set {/ a °° jp < oo}, the continuous local 

martingale f* ^j^ 3 - converges as t —>• oo. But we have seen that with unit probability In Rt 
does not converge. ■ 


We shall also use an absolute continuity result ([28], Exercise XI. 1.22 or [18], Proposition 
6.1.5.1). On the canonical space = C([0,oo), [0, oo)), we denote by R the canonical map 
R-tfuj) = u(t), by 7 Zt = cr(R s , s < t) the canonical filtration and by P^ the law of the Bessel 
process of dimension d >2 starting at r > 0. Then, 


( 8 ) 


P“ 


n t ~ 



(d - 2) 2 [* ds\ 2 

8 Jo R 2 s ) ' r 


Kt ■ 


3. Four useful tools 

The following simple lemma is a comparison lemma. It will play an ubiquitous role in our 
proofs. 

Lemma 4. For T > 0, a > 0, let x\ and x 2 be nonnegative continuous solutions on [0, T] to 
the equations 

xi(t) = ^W + a/o‘4) 
x 2 {t) = + <* ft Ztfj 

where v\, v 2 are continuous functions such that 0 < ui(0) < U 2 ( 0 ), and v 2 — v\ is nondecreas¬ 
ing. Then x\ (t) < x 2 {t) on [0,T]. 

Proof. Assume there exists t £ (0,T] such that x 2 {t) < x\ (t). Set 

t := maxjs < t : xi(s) < x 2 (s)} . 

Then, 

x 2 (t)-xi(t) = x 2 (t) -XI (t) + (v 2 (t) -V!(t)) - (v 2 (t) -Ui(r)) + a f*(^ _ ^)ds 

> 0 , 


a contradiction. 
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The following elementary lemma will also be repeatedly used. 

Lemma 5. Let Q(x,y ) = ax 2 + bxy + cy 2 be a second degree homogeneous polynomial. Then 
Q is nonnegative on the whole S if and only if a > 0, c > 0 and b > —2y/ac. 

Proof. 

Taking x = 0 (resp. y = 0) we see that a (resp. c) must be nonnegative. Then we rewrite 

Q(x, y) = ( y/ax - \[cy) 2 + (b + 2 ^/ac)xy. 

So Q is nonnegative for all x > 0 and y > 0 iff b > —2y/ac. ■ 


Another main tool will be the following convergence result, whose argument goes back to 
McKean ([211. p.31 and p.47). The statement and the proof below are borrowed from ([28]. 
Theorem V.1.7 and Proposition V.1.8). The only change lies in the introduction of a stopping 
time r up to which the local martingale M is now defined. 


Lemma 6. Let M be a continuous local martingale defined on [0, r) where t is a stopping 
time. Let ({M) t )o<t< T be its quadratic variation and let (M) T := lim < oo. 

(1) On {{M) r < oo}, Mt exists a.s. in M. 

(2) On {(M) t = oo}, limsup t _^ T M t = — lim inf t _>. T Mt = +oo a.s. 


Proof. (l)For any p > 1 let 

(j p = inf{f > 0 : (M)t > p}. 

In order every term in the following to be well defined we introduce a nondecreasing se¬ 
quence of stopping times (r n ) n >i with limit r such that each stopped process (M TnAt )t> o is 
a uniformly integrable martingale for any n > 1. The stopped process (M iATnAa - p )i >q is a 
L 2 -bounded martingale and as m, n —>• oo 

E [{M Tn A(T p - M T m Acr p ) 2 ] = ^[\{ M )r n Aap ~ (M) TmA(Tp \] ->• 0. 

We set 

m (p) ■- lim M t A(J = lim M tAa 

n —>oo t—*r 

On {(M) r < oo}, the stopping times a p are a.s. infinite from some p on and we can set 
M t := lim p _ KX) . Thus on this set 

.1 f —y AL t as t —y t . 


(2) Let for any t > 0 

T t = inf{0 < s < t : (M) s > t}. 

There exist an enlargement (0, (J}),P) of (f2, (J 7 ^),P) and a Brownian motion /? on Ll inde¬ 
pendent of M such that the process 


r M Tt if t< (M) t 

t \ M r + k-{M) r if t > (M) t 
is a standard linear Brownian motion. As 

lim sup Bt = — lim inf Bt = +oo a.s. 

t.—too t.—too 


we obtain on {(M) T = oo} 


lim sup Mr; = — lim inf Mr, = +oo 

t-y oo 


lim sup M t 

t—yr 


— lim inf Mt = +oo 

t—fr 


and therefore 
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We will also need the following consequence of the results in [2] on multivalued stochastic 
differential systems, completed with the method used in [3] and developed in m to check 
the lack of additional boundary process. 

Proposition 7. Let a > 0, 6 > 0, a = (cr*;z,j = 1,2) a 2 x 2-matrix, ( B,C ) a Brownian 

motion in the plane, b\ and b 2 two Lipschitz functions on M 2 , Zq and Z q two J-Q-measurable 

nonnegative random variables. There exists a unique strong solution (Z^Z 2 ) to the system 

Z\ = Zq + a\B t + cr\Ct + a Jo ;fj + Jo b\(Z\, Z^)ds 

z t = z o + °\Bt + cr%C t + 6 fg + f* b 2 (Zl , Zl)ds 

with the conditions Zj > 0 if 5 = 0 and Z[ > 0, Z 2 > 0 if 6 > 0. 


It is worth noticing that the solutions to © enjoy the Brownian scaling property. It 
means that if ( X , Y ) is a solution to © starting from (Xo, Yq) with driving Brownian motion 
( Bt,Ct ), then for any c > 0 the process (X' t := c~ 1 X c 2 t ,Yf := c~ 1 Y c 2 t ;t > 0) is a solution to 
© starting from (c 1 Xq,c 1 1q) with driving Brownian motion (c 1 B c 2 t ,c 1 C c 2 t ). 


4. Avoiding the corner 

We shall see in Section 5 that existence and uniqueness of the solution to © are easily 
obtained as soon as the solution process keeps away from the corner. Thus the question of 
attaining the corner in finite time is of great interest. For some class of reflection matrices, a 
necessary and sufficient condition is given in [29]. Unfortunately, we are not able to provide 
such a complete answer and we have to be content with a collection of sufficient conditions 
ensuring nonattainability. We will just see at the end of Section 5 a degenerated case where 
the corner is reached in finite time with full probability. 

Our sufficient conditions are stated in the following theorem. In some cases (conditions 
C 2a and C 2 b below), the comparison with Bessel processes suffices to conclude. In other cases 
(conditions C\ and C 3 below), we are first looking for a C' 2 -function / on the punctured 
quadrant S° with limit —00 at the corner. Then we use Lemma E] to show that f(X t ,Y t ) 
cannot converge to —00 in finite time. 

Theorem 8. Let (X,Y) be a solution to ©). We set 

r° := inf {t > 0 : (X u Y t ) = 0 } 

with the usual convention inf 0 = 00 . Then P(r° < 00 ) = 0 if one of the following conditions 
is satisfied: 

(1) C\ : f3 > 0, 7 > 0 and — 1 < p < a + 5. 

(2) C 2a : a > \ and (3 > 0 

(3) C 2 j, : 6 > \ and 7 > 0 

(4) C 3 : There exist A > 0 and p > 0 such that 

• Xa + p'-j > 0 

• A/3 + p5 > 0 

• (\/A(Act + p 7 ) + \Jp{ A/3 + pd))~ > + p~ + 2p\p). 

Proof. 

Condition C\. For e > 0 let 

= 1 {(ATo,y 0 )=o} inf {t > 0 : X t + Y t > e} 
r°’ e = M{t > a £ : (X t ,Y t ) = 0 } . 




8 


A TWO-DIMENSIONAL OBLIQUE EXTENSION OF BESSEL PROCESSES 


As e | 0, a e | 0 and r°’ e } r°. We set Rt = Af 2 + Y 2 . From Ito’s formula we get for 
t €E [cr e , r 0,e ) 

In Rt 

= In R„ f + 2 £ XdB s +Y s dC s + 2 £ + 2 Jl (/?§ + 7 £)£ - 4 p /* e 

> In i? CT e + Mi + 2 —^jp—^ds 

where M is a continuous local martingale on [(x e ,T 0,e ) and P\(x,y) is the second degree 
homogeneous polynomial 

Pi (x, y) = {a + < 5 ) ( x 2 + y 2 ) - 2pxy. 

Using Lemma [5] we check that this polynomial is nonnegative on S if p < a + 5. Therefore 

P i (X s , Y s ) 


0 < 


R 2 


-ds < oo 


From Lemma [ 6 ] we know that depending on whether (M) r o, e is finite or not, the local mar¬ 
tingale Mt either converges in 1 as t -> r 0,<E or oscillates between +oo and —oo. It cannot 
converge to —oo. Thus R T o, e > 0 on {r 0,e < oo}, which contradicts the definition of the 
moment r 0,e . This proves that r 0,<E = oo for every e > 0, and therefore, r° = oo a.s. 
Condition C 2 a (resp. C-jb)- Recall the Bessel processes U and V in ([2]). From Lemma 1 we 
get X t > Ut (resp. Y t > Vt) and in this case 0 is polar for U (resp. V), so X t > 0 (resp. 
Y t > 0) for t > 0. 

Condition C$. We use again the notation for e > 0, a e and r 0,e . We set St = XX t + yY t for 
t > 0, A > 0 and y > 0. From Ito’s formula we get for t € [a £ , t 0,£ ) 

In S, 

= In S a e + XdB sM d Cs + (Aa + yy) f* c + (A/3 + yd) f* c 
-\{X 2 + y 2 + 2p\y) || 

= In So-+M t + £%$$& ‘ 

where M is a continuous local martingale on [(T e ,r 0,e ) and P 2 (x,y ) is the second degree 
homogeneous polynomial 

P'iixyy) 

= A (A/3 + yd)x 2 + y(Xa + yy)y 2 + [A(Aa + yy) + y(X/3 + yd) — ^(A 2 + y 2 + 2pXy))xy . 
Using again Lemma 0 we see that P 2 is nonnegative on S if 


-(A 2 + y 2 + 2pXy) — [A(Ack + yy) + y(X(3 + yd )] < 2 \J A(A/3 + yd)y(Xa + yy)- 


This is exactly the condition C3. Therefore 


0 < 


/' 


P2(X s ,Y s ] 

X s Y s S 2 


■ds < 00 


and so 


0 < 


/ 

J cr 


-0,e 


P2(X s ,Y s ) 

XsYsS* 


ds < 00 . 


Similarly, using Lemma [ 6 ] again, we see that the continuous local martingale M either con¬ 
verges to a finite limit or oscillates between +00 and —00 when t —> r 0,e . It cannot converge 
to —00 and thus S T o, £ > 0 on {r 0,e < 00 }, proving P(r 0,e < 00 ) = 0. Letting finally e —> 0 we 
obtain P(r° < 00 ) = 0 . ■ 


Condition C 3 is not explicit. We give two concrete examples when this condition holds 
true. 














A TWO-DIMENSIONAL OBLIQUE EXTENSION OF BESSEL PROCESSES 


9 


Corollary 9. Assume p = 0 
( 10 ) * 


a = 5 and |/3| = | 7 |. Then the condition C 3 is satisfied if 

ft 2 < a — | when (3 = —7 
—/3 < a — ^ when ft = 7 < 0 . 


Proof. In both cases we take A = p > 0. 
When f3 = —7 > 0, the condition C 3 writes 


• a — fi > 0 

• 1 < (V« + /3 + — /3) 2 = 2a + 2-^/a 2 - /3 2 . 

If /3 2 < a — |, then 

a 2 — /3 2 > a 2 — a + 7 > 0 
4 - 

and 

(1 - 2a ) 2 = 1 - 4a + 4a 2 < 4(a 2 - /3 2 ). 

When fi = 7 < 0, the condition C 3 writes 

• a + fi > 0 

• 1 < (2-^a + /3 ) 2 = 4(a + fi) 

and this is —fi < a — ■ 


Corollary 10. Assume 
• max{a, 5} > \ 

Then P(r° < 00 ) = 0 . 

then condition C* 2 a holds true and the conclusion 
We take A = 5 and p = —ft. Then A/3 + p5 = 0, 

da — fi'y 

da — /3a(2p — ^) 
f(5 2 + /3 2 - 2p(3d) 

W~PP ) 2 

0 


Proof. We may assume a>\. If /3 > 0, 
follows. If /3 < 0, we will use condition C 3 . 

\a + p'y = 
> 

> 

> 


and 


(\/ A(Aa + /j. 7 ) + \J //(A/3 + pd)) 2 — |(A 2 + p 2 + 2p\p) 
= \{\a + p'y) — \(X 2 + p 2 + 2p\p) 

> a(5 2 +/3 2 -2p/M)-l(h 2 + ) 0 2 -2 / 9/3J) 

> (a — ^)(<5 — p(3) 2 

> 0 . 


5. Existence and uniqueness 

We now proceed to the question of existence and uniqueness of a global solution to (H|). 
We consider separately the three cases: first /3 > 0 and 7 > 0, second fi'y < 0, third /3 < 0 
and 7 < 0. In the first and second cases, we construct the solution by switching from one 
edge to the other and patching the paths together. Thus it is essential to avoid the corner, 
as it was supposed in |30j in order to weakly approximate an obliquely reflected Brownian 
motion. The third case uses a different proof and does not requires avoiding the corner. All 
three proofs heavily use the comparison method of Lemma 01 
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5.1. Case j3 > 0 and 7 > 0. 

Theorem 11. Assume fi > 0, 7 > 0 and one of the conditions C\. C^o, C 2 b, C3 is satisfied. 

( 1 ) There is a unique solution to HP in S°. 

(2) There is a solution to HP in S starting from 0. 

(3) If aS > fi'y, there is a unique solution to HP in S. 

Proof. 1. Let a > 0, e > 0 and define for (x, z) € M + x R 

if e (x,z):= 


( 11 ) 


max (72; + z, ae) 

This is a Lipschitz function. From Proposition [7] we know that the system 
Xf = X 0 + B t + a + a [3 Zf)ds > 0 

Z\ = -7X0 + a{Yg + l{y 0 =o}o) - l B t + aC t + a(ad - fi'y) ^fi e {X e s ,Z e s )ds 
has a unique solution. Let 

Ty ■= inf{t > 0 : + Z\ < ae} . 

IfO< 77 <e<awe deduce from the uniqueness that ( X e , Z e ) and (X' n , Z v ) are identical on 
[0, Ty]. Patching together we can set 

X t := lim^o Xf 

Yt ■= Um e ^ 0 i(7 Xt + Zf) 

on {(oj,t) e!lx [0, 00 ) : Yq(uj) > 0 and 0 < t < Ty(w)}, where 

Tv := lim tL . 

e—10 

On this set, (X,Y) is the unique solution to (JTJ). As we already noted, we have X t > Ut and 
Yt. > Vt . Therefore, on {Yo > 0 } n {Ty < 00 }, 


fo° & < JtT" * < 00 and jf * < # < 00 


and we can define 

( 12 ) 


X 0 

r Y 


~ = X 0 + B r o + a ff % +13 ff 


0 Y. 

Yr° ■■= = Yq + Cj +Tfa° 3 T+SJo° $;■ 


‘Y L 7 ' Y ‘Y 

We have Y r 0 = 0. From Theorem Owe know that X T o > 0. In exactly the same way we can 
construct a solution on {Yo > 0} in the interval [Ti,T 2 ], where T\ = Ty, T 2 = infjt > T\ : 
Xt = 0}. Iterating, we get a solution on {Yo > 0} x [0, lirrq^oo T n ) where 

Tip := inf{t > T 2p -\ : X t = 0} 

T 2p + 1 ■= inf{t > T 2p : Y t = 0} . 

On {Yo > 0} O {Hindoo T n < 00 } we set X^^^ := Hindoo X T2p = 0 and Y^^ Tri := 
liirip^oo Yr 2p+1 = 0. The polarity of 0 entails this is not possible in finite time and thus 
lim n _ ) . 00 T n = 00 . So we have obtained a unique global solution on {Yo > 0}. In the same 
way we obtain a unique global solution on {Xg > 0} and as P((X 0 , Yg) = 0 ) = 0 the proof is 
complete. 

2. Assume now Xg = Yg = 0. Let (y n )n> 1 be a sequence of real numbers (strictly) 
decreasing to 0. From the above paragraph it follows there exists for any n > la unique 
solution (X n ,Y n ) with values in S° to the system 

Xf = B l + af‘$ + l3f‘$ 


yn 


— Un + C t + 7 fg jfTT + d /(- 


ft ds_ 

0 
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Let 

r := inf {t > 0 : X t n+1 < X?} . 

Using Lemma [4] we obtain Y t n+l < Y t n on [0 ,t]. We note that (X", Yfi) £ S° on {r < oo}. 
On {Yfi +1 = y T n }n{r < oo}, since X™ +1 = X" and the solution starting at time r is unique, 
it follows that X” +1 = X™ and Yfi +1 = Y t n on [r, oo). On {Yfi +l < Yfi} n {r < oo}, the 
continuity of solutions at time r entails there exists p > 0 such that Yfi +1 < Y t n on [r, t + p]. 
A second application of Lemma 0] proves that X”" 1 " 1 > X" on [r, r + p], a contradiction to 
the definition of r. Therefore P(r = oo) = 1. It follows that X t n+1 > Xf and Y t n+1 < Y f n for 
any t £ [0,oo), and we may define 

X t := lim t X } 1 Y t := lirn | Y t n . 

n —^oo n—>oo 

As Y t n > Vt where (U, V ) is the solution to ([2]) with Xq = Yq = 0, we have 
Xt = B t + a lim n _^ 00 f* + fi lim^oo yk 

= B t + a f 0 + ft Jo 

< oo 

and also 

Yt — lim n ^ 00 y n + Ct -\- liniT^—^oo Jq + (5 linin—^oo Jq y^r 

= C« + 7/„‘^ + 4/o‘# 

< OO . 

3. Assume finally a<5 — /Ly > 0 . As the conclusion holds true if (3 = 7 = 0 , we may also 
assume f3 > 0. Let (X,Y) be the solution to (0|) with Xo = Yo = 0 obtained in the previous 
paragraph and let (X',Y') be another solution. Considering ( X n ,Y n ) again and replacing 
(X n+ \Y n+1 ) with (X',Y'), the previous proof works and we finally obtain X[ > Xt and 
Y( < Y t . Then, 

(13) 0 < 5(X[ - X t ) - (3{Y{ ~ Yt) = f\aS - /Ly)(^ “ Y~) ds < °- 

Thus X[ = X t and Yf = Y t , proving uniqueness. Replacing (5, (3) with ( 7 , 0 ) in equation (fTUl) 
we obtain the same conclusion if 7 > 0 . ■ 

Remark 12. The statement in Theorem I 111 is not complete since the problem of uniqueness 
when starting at the corner and ad < /Ty is not solved. When considering the solution (X, Y) 
in the above proof of existence, we have noted that X’ > X and Y' < Y for any other 
solution (X',Y'). Thus uniqueness in law would be sufficient to obtain path uniqueness. A 
possible way to prove weak uniqueness could be the method in pLj. This would be far from 
our trajectorial methods and we don’t go further in that direction. 


5.2. Case /Ly < 0. 


Theorem 13. Assume /Ty < 0 and one of the conditions C^a or C 3 is satisfied. Then, there 
exists a unique solution to m in s°. 

Proof. Assume first (3 > 0, 7 < 0. The proof is similar to the proof of 1 in Theorem 1111 
The only change is that now Y t < Vt . Therefore, on {Lo > 0} n {fiy < 00 }, 


. [X ds [X ds 

L YXx- Y °- c> r-y 0 w. 


< 00 
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and we can define X T o and Y T o as previously done. The application of Theorem [H] to the 
process on the time interval [0, Ty] shows that X T o > 0 and we can iterate the construction 
as in Theorem I 111 The proof if (3 < 0, 7 > 0 is analogous. ■ 


5.3. Case (3 < 0 and 7 < 0. In this case we can give a full answer to the question of existence 
and uniqueness. When \p\ < 1, our condition of existence is analogous to the condition found 
in m for the reflected Brownian in a wedge being a semimartingale, i.e. there is a convex 
combination of the directions of reflection that points into the wedge from the corner. It 
amounts to saying that the interaction matrix 

a (3 
7 5 

is completely-5 in the terminology of ESI El US EES ESS EE- 

Theorem 14. Assume (3 < 0 and 7 < 0. 

(1) If ad > f3'y, there exists a unique solution to f7]l in S. 

(2) If ad < /?7 and 1 + p + \a + 7 I + |/3 + <5| >0, there is no solution. 

(3) If l + p = a + r y = /3 + 5 = 0 and (Xq, To) / 0 there exists a unique solution. 

(4) Ifl+p = a + 'y = /3 + d = 0 and (Xq. Yq) = 0 there is no solution. 


Proof. 1. Assume first ad > /3 r y. 

a) Existence. Let ( h ni n > 1) be a (strictly) increasing sequence of bounded positive nonin¬ 
creasing Lipschitz functions converging to 1/x on (0, 00 ) and to +00 on (— 00 ,0]. For instance 
we can take 

h n (x) = (l-^)y on [i,oo) 

= n - 1 on (- 00 , i]. 

We consider for each n > 1 the system 

(u) W = X 0 + B t + af*-$k + /3fih n (Y?)ds 

Y t n = Fo + Cf + 7 Jq h n (Xg )ds + d Jq ^ . 

From Proposition [7] it follows there exists a unique solution to this system. We set 

r := inf{s > 0 : X(f +1 > X(f}. 

We have h n+ \ (Xf +l ) > h n (Xf) on [0, r]. A first application of Lemma 0] shows that Y t n+1 < 
Y t n on [0, t\. Since h n+ i(Yf- +1 ) > h n (Yf-) on {r < 00 }, we deduce from the continuity of 
solutions that there exists a > 0 such that h n+ i(Yf l+1 ) > h n (Y t n ) on [r,r + a]. A second 
application of Lemma 0] shows that Xf +1 < Xf on [r, r + a\, a contradiction to the definition 
of r. Thus P(r = 00 ) = 1 proving that on the whole [0, 00 ) we have Xf +1 < Xf and 
y-n+l ^ yri' 'pj len W g can gg^ f Qr an y £ g [ 0 , 00 ) 

X t := lirn X? and Y t := lim Y t n . 

n^-oo n—> OO 

If ad > P'y, there is a convex combination of the directions of repulsion pointing into the 
positive quadrant, i.e. there exist A > 0 and p > 0 such that Xa + p'y > 0 and pd + Xf3 > 0. 
For n > 1 and t > 0, 


(15) 


XU t + pV t > A X? + pY t n 

> AA 0 + pYo + A Bt + pC t + (Aa + p'y) f 0 + (pd + A/3) f 0 -yr ■ 


Letting n —>• 00 in (fl5l) we obtain 


A TWO-DIMENSIONAL OBLIQUE EXTENSION OF BESSEL PROCESSES 


13 


Then we may let n go to oo in (TH1) proving that (X, T) is a solution to CD), 
b) Uniqueness. Let be another solution to ([[]). Replacing (X n+1 ,Y n+1 ) with (X',X') 

we follow the above proof to obtain for t £ [0, oo) and n > 1 

Xl<X t n and Y/<Y t n 

Letting n —>• oo we conclude 

X[ < X, and Yj < Y, . 

With the same A > 0 and p > 0 as above, 

0 < A(X t - X') + p(Y t - YD = f\(Xa + --]-) + (p8 + A/3)d - ^)}ds < 0 

JO r s * s 

and therefore X[ = X t , Y{ = Y t . 

2. If aS < fX/ there exist A > 0 and p > 0 such that Act + p'y < 0 and p5 + A/3 < 0. For that, 
just take 

( 16 ) “ '“'-V 


a u 
— < L < 

—7 A 5 


Let us consider the nonnegative quadratic (A 2 + / u 2 + 2pA / u). It is positive if p > 0. It is larger 
than (1 — p 2 ) A 2 > 0 if — 1 < p < 1. It may be positive if p = — 1 and \a + y| + \fi + <5| >0 
because from (fl 6 l) we may take A 7 ^ p. Thus, if (X,Y) is a solution to flTJ), for any t> 0, 

0 £ XX{ + pY^ ft XXq + P^o T A Bt + pUt. 

This is not possible since the paths of the Brownian martingale {\ 2 +p 2 +2p\p )~ 1 / 2 (\B t +pCt) 
are not bounded below. So there is no global solution. 

3. Assume now 1 + p = \a + q| + \/3 + 5| =0 and Xo + Yq > 0 . The system becomes 

X, = x 0 + B t + aj‘f ; -Sf‘^ > 0 

Y, = + > 0 . 

This entails for any t > 0 

X t + Y t = X 0 + Y 0 

and the first equation in (tl7]) reduces to 


(17) 


(18) 


0 < Xt — Xn + Bf\- ot 




ds 


X 0 + To - X s 


< X 0 + T 0 , 


Clearly pathwise uniqueness holds for equation (1181) since if there are two solutions X and 
X for any t > 0 

(X t — X ' t ) 2 = 2 a f 0 (X s — X')(^T — Yr)ds — 26 f Q (X s — X')( x 0 +Vo-X s ~ X 0 +Yb-X' s 

0 . 


Consider now for any c > 0 and 0 < cq < c the equation 


(19) 


0 < Z t = c 0 + B t + a 




ds 


< c. 


C - Z s 

The solution Z is a Brownian motion perturbed by a drift deriving from the concave potential 

alnz + 6 ln(c — z) for 0 < z < c. 

From |2j we know that equation (1191) has a unique solution living on the interval [0, c]. Thus 
there is a weak and then a unique strong solution to (fT 8 l) . Setting 


Y t = X 0 + T 0 - X t 
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we obtain a unique strong solution to (HD- 
4. We must have for any t > 0 

X t + Y t = X 0 + Yb = 0 

and thus X t = Y t = 0. But this is not possible for a solution to (fTj) . ■ 

We end this section by considering a degenerate case where the solution hits the corner 
with probability one. 

Proposition 15. Assume p = 1, ad > /Ty, max{a;,<5} < I and max{/3, 7 } < 0. Then, 
P(t° < 00) = 1 . 

Proof. We set 

Ty := inf{i > 0 : X t = 0} 

Ty := infjf > 0 : Y t = 0} . 

The dimension of each Bessel process U and V in ([2]) is less than 2, and X < U, Y < V. 
Then P(r^ < 00 = 1) and P(ry- < 00 = 1). Assume first a + /3 < 7 + 5. Using Lemma 0] we 
obtain Xt < Yt < Vt on the time interval [t^cxd). Therefore P(r° < 00 ) = 1. Same result 
when a + /3 > 7 + <5. ■ 

In the following pictures, we display the directions of interaction r x and r y defined in (0|) 
in three illustrative instances. 




6. Avoiding the edges 

We now consider the question of hitting an edge of the quadrant. Remember the definitions 

. := inf{t > 0 : X t = 0} 

{ ’ Ty := inf{t > 0 : Y t = 0} . 

We already know that P(ty <oo)=0ifa>^ and (5 > 0. Conversely, a comparison with 
the Bessel process U shows that P(r^ <oo) = lifa<^ and /3 < 0. If we know that the 
corner is not hit and a > we can get rid of the nonnegativity assumption on f3. Since we 
are only interested in one coordinate, we are looking for a function that is C 2 on ( 0 , 00 ) and 
goes to —00 when approaching 0. A natural candidate is the logarithmic function. 

Proposition 16. Assume P(r° < 00 ) = 0. If a>\, then P(ty < 00 ) = 0. 

Proof. For 7 > 0 let 

6 x = 1 {X 0 = 0 } inf{t > 0 : x t >rj} 

Ty v = inf{i > 6 \ : X t = 0} . 
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As r) | 0, 0\ 0 and 4- t y- F° r ^ € [^x’ r x ?7 )> from Ito’s formula 

( 21 ) lnX t = lnX g ^ + ^ + (a - 5 ) jf| 

Since P(r° < 00) = 0 , on the set {t^ v < 00} we have Y o, v > 0 because X r o, v = 0 . On this 
set, from the definition of an 02 BP, 


P fe^ x s Y a ■ 


0,77 

r* ds_ 
Xs 

0,77i 


< 00 


0,77 

[ T X ds 

/ v < °°' 
lei Ys 


As Y s > 0 on some interval [x, T x\ with positive measure and X s > 0 on [ 9 V X , x], we see that 


m 


0 ^ 
°x 


ds 

XsY 


< 00. 


As t —>• Ty* 1 , the local martingale in the r.h.s. of (12111 cannot converge to —00. This entails 
that P (rff 1 < 00) = 0 and therefore P(r^ < 00) = 0 . ■ 


We may also be interested in hitting either edge. This time we are looking for a function 
that is C 2 in the interior of S and goes to —00 when approaching either edge. 

Proposition 17. Assume a > \ and S > Then P(r^ < 00 ) = P(ry- < 00 ) = 0 if one of 
the following conditions is satisfied: 

( 1 ) /3 > 0 

(2) 7 > 0 

(3) 0</3y< (a - i)(<y-i). 

Proof. The proofs in the cases (1) and (2) ) are direct consequences of Theorem [ 8 ] (under 
conditions C^a or C 26 ) and Proposition [T 6 j Assume now the conditions in (3) hold true. For 
e > 0 let 


<7 e = 1 


=o> inf{t > 0 : X t Y t > e] 


{X 0 Y 0 = 0} 
r e = inf(t > cr e : X t Y t = 0} . 

For A > 0 and fi > 0 we set 

Rt = A In X t + ft In Y t . 

From Ito’s formula we get for t € [cr e ,r e ) 

R t = Ra‘ + fie( yi- dB s + frdCs ) + fie[ A + ^ X xIys^^ s 

= Rat+Nt + fle S yf } ds 

where N is a continuous local martingale and P^(x,y) is the second degree homogeneous 
polynomial 

P-i(x, y ) = ~ ]jf)x 2 + A (a - i )y 2 + (A/3 + ni)xy. 

If the conditions (3) are satisfied, we may take 

A = 2(a — f)(S — \) — /3y > 0 
h = P 2 

Then, using again conditions (3), we check that 
Pz{x,y) = /3 2 ( 5 - ^ x 2 + 


2 a- 


1 


5 - \ ) - h 


> 0 . 


1 


a--)y +2 [a- 




is nonnegative on S. The proof terminates as previously in Theorem [ 8 ] . 
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7. Associated local martingales 


We easily check that if max {/3, 7 , aS — /5^y} > 0 and if p > — 1, there exist some A > 0 and 
p > 0 with A + > 0 such that XX t + pY t is not less than A(Xo + B t ) + p(Yq + Cf) which is 

proportional to a real driftless Brownian motion. Therefore 

limsup(AA t + fiYt) = + 00 . 

i—>• 00 

An usual way in the study the asymptotic behavior of the solutions to stochastic differ¬ 
ential equations is introducing associated martingales. This is carried out through scale 
functions ([28] Section VII.3). There is no equivalent functions on the plane. However, in 
some particular cases we can find simple functions of 02BPs that are supermartingales or 
local martingales. 


Proposition 18. Assume the following set of conditions: 



a > 1/2 

(22) 

5 > 1/2 

P + 7 > 0 

25-1 1 2a—1 — P 


p > -1. 

Then 

M t := x]- 2a Y^- 2S 


is a positive supermartingale on (0, 00 ) which tends to 0 as t —>• 00 . It is a local martingale 
if the third inequality in |jjj| ) is an equality. 


Proof. For e > 0 let 

<7 e = l{x 0 v 0 =o} inf {^ > 0 : x t Y t > e} 
r e = inf{7 > <r e : X t Y t = 0} . 


As e | 0, a e I 0. Up to r e , we get from Ito’s formula on {a e < 00 } 


j\/f j\/f 1 (1— 2<y.)Y s dB s -\-( < l—25)X s dC s 

M t ~ + J a e - X ‘2 a y‘2S - 

+[(1 - 2 a)/3 + (1 - 26h + p(2a - 1)(25 - 1 )] ^0^. 

The function f(x, y) = x l ~ 2a y l ~ 2S is C 2 in the interior of S and goes to +00 when approaching 
the edges of the quadrant, whereas the finite variation part in the semimartingale decompo¬ 
sition of Mt is nonincreasing. We apply again Lemma [ 6 ] and obtain that P(T e < 00 ) = 0. 
Then, letting e —>• 0, we see that M is a positive supermartingale on (0, 00 ). As such it tends 
to a limit H > 0 when t —>• 00 (|28j. Corollary II.2.11 ). If /3 > 0 and 7 > 0, then X > U 
and Y > V where U and V are the Bessel processes in ([2]). We have seen in Section 2 that 
Ut —» 00 and Vt —^ 00 as t — y 00 ; so M t —>• 0 as t —>• 00 . If now min{/ 3 , 7 } < 0, say 7 < 0, 
then we have Y < V. From Corollary [3] we deduce that for any e > 0 


(23) 



We consider the quadratic variation (M) £ given by 


(M)t 



[(1 - 2a) 2 Y7 + (1 - 2 d) 2 Xf + 2p(l - 2a)(l - 26)X s Y s ]ds 


+ 


(1-2 <5) 2 

YT~ 


+ 2p 


(1—2a)(l—25) 
X.Y. 


]ds. 



M 2 


(1 - 2 5) 2 
Y 2 

S 


ds 


If p > 0, 
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and if p 2 < 1 , 

(M)t > (1 - p) £ 

In both cases, using rt23|) we should have (M)^ = oo on the set {H > 0} n {a e < oo}, and 
then limsupj M t = — liminf t M t = oo, a contradiction with M t —>• H. Thus H = 0 a.s. ■ 

During the proof we have seen that P(ry < oo) = P(ry < oo) = 0. In fact this is not a 
new result since here the conditions in Proposition [17] are in force. There is another case of 
interest. 

Proposition 19. Assume a = 5 = 1 /2, pf3 7 < \P'y\, and there exists a solution to (Q]) 
satisfying P(r° < 00 ) = 0. Then 

M t :='ylnX t -plnY t 

is a continuous local martingale on (0, 00 ) and 

lim sup Mt = — lim inf Alt = 00 

t— KX> t—>OQ 


Proof. We deduce from Proposition [T3 that X t > 0 and Y t > 0 for any t > 0. Let again 
for e > 0 

= l{x 0 y 0 =o} inf {t > 0 : X t Y t > e}. 

Now Ito’s formula gives 

rt ( dB s dC s 

X s P Y s 


Mt = M a e+ I ( 7 -^r- - P 
and the quadratic variation is 

(M)t = [ 

J a 

Then, if p/3 7 < 0, 


r ' ( 7 2 , p 2 0 pt ,, 
U .; 2 n 2 p x s y s Us - 


(M) e t > [ 

J (J 




and if p 2 < 1 , 


{M)t > (1 - P 2 ) max | J * j^ds, j ^ . 

If P > 0, then X > U where U is the Bessel process of dimension two in ([2]) and 


If P < 0, then X < U and 


lim sup (In X t ) > lim sup (In Ut) = 00 

> 00 t—yoo 


lim inf (In Xf ) < lim inf (In Ut ) = —00 

t—> OO t—>00 


But 


and 


In X t = In X a e + 




r* ds 1 r 
jXsYs - 2 J a AX 2 ' Y? 


ds 

a e Yv s J a e X S Y S 

t 


1 1 . , 

T 7 T 7 T ) ds. 


(T fc \-“-s S 


Thus In X f would converge a.s. as t —>• 00 if 

p0 ° ds 


f 

J (J 


X 2 

(7 C - /Y S 


< 00 and 


f°° ds 

L Y F 


< 00. 
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Therefore (M)^ 


oo on {a e < 00 } and the assertion is proved. 


8. Absolute continuity properties 

In this section we suppose p = 0. In some cases we easily obtain an absolute continuity 
property between the laws of 02BPs with various parameters. When there exists a unique 
solution to (P), we denote by Ff’y’ 7,S the law on C'(R+, S) = C'(M + ,M + ) x C(M+,M+) of the 
solution starting at (x,y) € S. We denote by (U,V) the canonical map (Ut(u,v),Vt(u,v)) = 
(u(t),v(t)) and by lit = &((U S , V s ), s < t) the canonical filtration. Recall that we denote by 
P(? the law of one-dimensional Bessel process of dimension d starting at r > 0. 

Proposition 20. Assume p = 0, 5 > | and y > 0. Then 

(24) KS” l».= -p {a f‘ y -<#f 0 4 k, - f J ‘#} •- p2 ° +l ' 18p “ +1 • 


Proof. Under P ^ Q+1 ® F 25+1 the process 


B t := U t - x - a 


r* ds_ 
Jo ~Us 


is a one-dimensional Brownian motion. The assumptions on 6 and on y imply that V) > 0 
for any t > 0 and thus 

ds p2<5+l 


v 2 (s) 


< 00 


— a.s. 


Therefore, for P ^ +1 -almost every v, is a P^ Q+1 -centered Gaussian variable with vari- 

ance Jo ipfe and 


exp < f3 


dB s 

v(s) 


P f 

2 Jo 


We see that 


Z t := exp \ /3 


ds 

v(s) 2 
* dB s 


dP: 


i 2 a+l 


W ) dF 25+ \v) = 1. 


V, 


fP ft ds_ 

^ Jo V? 


is a P ^“ +1 ® P^ +1 -positive martingale with expectation 1 . Setting for any T > 0 


Q t :=Z T .F 2 x a+1 ®F 2 y 5+1 \ Ut 


we check that under 


Ut — x — a 


Lw.< 


ds 

V' 


0 <t<T 


is a real Brownian motion independent of the Brownian motion {Vt — y — 6 Jq y, 0 < t < T}. 
Thus, 

Qt = Kf 0 - 6 K ■ ■ 


A second set of conditions is obtained using Novikov’s criterion. 
Proposition 21. Assume p = 0 and the following set of conditions is satisfied: 

\P\ <S~ 1/2 M<«-l/2 

x > 0 y > 0. 


( 25 ) 
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[* ds 
'o UsVs 


The process 

Z t := exp | (3 f ^-7 f ^ ~ W + 78 ) 

l Jo v s Jo Us Jo v s v s z j 0 Vg 

is a P^“ +1 <8) P 2 y S+1 -positive martingale with expectation 1 and 
(26) F^’ S \ Ut =Z t .F 2 x a+ 1 ®F 2 y 5+1 \ Ut . 

Proof. Under P^ Q+1 <g) P^ +1 the processes 


^ f t ds__l ^ f* ds_ 

T J 0 v?~Y J 0 u! 


R rr f ds 

B t ■= U t - x - a — 

.In Us 


and 


C t := V t — y — 5 


L 


ds 
0 K 


are independent Brownian motions. Then Zt may be written 

r-t^o d An p2 rt ds 

^ Jo V? 

which shows that Zt is a positive local martingale. Using (j8]) we compute 


., 1 dB p dC, 

z, = exp ^/„ W + U„ 07 


^ f* ds_ 


])2a+l 

a; 


//exp{y/ 0 '|i}(iPj“ +1 ®(iP!f +1 = / exp {y dP; 

< OO 

since a Bessel process of dimension two has finite moments of any order. Finally, 

/ / «P {f /J f} exp $} dP? +1 « dPf +1 

= (/ exp {f f‘ f} <ff? +1 )(/ exp {£ /„* f } <JPf +1 ) 

< OO 

and Novikov’s criterion (ESI, Proposition VIII. 1.15) proves that Zt has expectation 1 with 
respect to P ^ Q+1 (giP ^ 1 . We easily see that {B t — f3 /q fr, 0 < t < T} and {Ct — 7 /)[' 0 < 

t < T} are independent Brownian motions under the probability with density and this 
proves that the new probability is P“ f 7, ‘ I U T - ■ 


9. Product form stationary distribution 
For a > 0 and c > 0, let T(a,c) be the probability measure on [0,oo) with density 


70W) := — x ^ e -™ 

r(a) 


and characteristic function 
(27) 


/■OO 

0(A;a, c):= / e lXx 'y(x-,a,c)d: 

Jo 


'x = 11 - 


i\ 


/ 0 V c 

For a > 0, 6 > 0, let 5(a, 6) be the probability measure on [0,1] with density 

P(x-a,b) := 

r(a)r(6) 
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It is well-known that the function g{x) = x d_1 is an invariant measure density on [0, oo) 
for the Bessel process of dimension d. To get a stationary probability we have to introduce 
a negative drift that entails positive recurrence. Then, the process 


X t 


Xq + B t + 



0t, 


with 6 > 0, has 7 ( 37 ; d, 20) as stationary density. As for the Brownian motion reflected at 
0 with constant drift —0, it has the stationary exponential density 2#exp{— 20x}. In the 
bidimensional case, the drifted obliquely reflected Brownian motion has a stationary density 
in the form of product of two exponential densities if and only if it satisfies at once DSHD]: 

• an invertibility condition on the reflection matrix; 

• a positivity condition on the exponential coefficients; 

• a skew-symmetry condition. 

Therefore it is natural to ask wether one can find similar conditions for drifted 02BPs ensuring 
existence of a stationary distribution in the form of product of two gamma distributions. 
The answer is positive and given in the next theorem. In fact, this may be considered as a 
consequence of the study in [23] which introduced a generalised reflected Brownian motion 
associated with a potential U regular on the whole real line. The difference is that here the 
logarithmic potential is defined only on the positive axis. The proof below is an adaptation 
of the proof in [23] to this special case. 

We introduce an additional constant drift (— 0 , — 7 ) in order to make the solution a recurrent 
process in the nonnegative quadrant. We consider the system 

m V, = x o + B l + af‘^ + /)f‘^-0t 

1 1 y, = r„ + c t + 7 / 0 , ^ + j/„‘g-»‘ 


with the conditions Xt > 0, Yt > 0. Changing probability through a Girsanov transformation, 
we easily check that Theorem [S] in Section 4 is still valid for this drifted system. Moreover, 
the proofs in Section 5 do not bother whether the Brownian motions B and C are drifted or 
not drifted. Therefore existence and uniqueness results in Section 5 hold true for the solution 
to ( 1281 ) . 


Theorem 22. Assume there exists a unique solution to l[28\) in S° or S. This process has 
an invariant distribution in the form T(a, c) < 8 ) T(b,d) with a > 1 and b > 1 if and only if at 
once 


• ad — /?7 7 ^ 0 

• arj — 70 > 0 


lP = i + i 


(invertibility of the interaction matrix) 
and 50 — /3r) > 0 (positivity of exponents) 

(skew-symmetry) 


Under these conditions, the unique solution is given by 


• a 

• b 


• c 

• d 


1 + 2 a 
1 + 25 


2a 

25- 


SO—(3rj 
aS-fi^y 

aif—^6 

aS—fi 7 * 


Proof. We first remark that under the above conditions a5 — (3 'y cannot be < 0, because 
in that case there is no solution if (3 and 7 are < 0 and if (3 and 7 are > 0 the skew-symmetry 
condition is not satisfied since 2p < 2 < (3/5 + 7 /a. So we could have written a5 — /Ty > 0 
as first condition. Let now 


q(x, y) = 7 (x; a, c) 7 (y; b,d) for x > 0 , y > 0 . 
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The infinitesimal generator of the diffusion 


1 d 2 d 2 <9 2 a B d 'yd d 

L = -(—2 + 2 + 2/9 ) + (-1- 6)— -h ( 1-??) R - 

o Qyz dxoy x y ox x y oy 


is given by 

d_ 
dx 


2 v dx 2 oy* oxoy' ' x 

Assume (Xq. Yq) has density q and characteristic function 

E [ e i(AX 0 + M Yo)] = Qj c ) ^ b, d ) . 

We set f(x,y ) = e l (*®+M2/). We want to prove that for any t> 0, A, and /r, 
E[f(X t ,Y t )] = E[f(X 0 ,Y 0 )} (= 0(A;o,c)0(/x;6,d) ) . 

It is enough to prove that 


| /*oo 

/ Lf(x,y)q(x,y)dxdy = 0 
2 o 


for any A and y. Let 


-R(A,//) := [-^(A 2 + /r 2 + 2pA/i) + iA(— + — — 0) + iy(— + - — r/)] 
2 x y x y 


and compute 

fo° fo° L f( x ’ y) dxd y = fo°° fo° R ( X ’ y)e^ Xx+ ^y{x\ O. c)y(y; b, d) dxdy 

= S(X,y)4>(X;a,c)4>(y,b,d) 

where S(X, y) is a second degree polynomial. In the last computation we used formula (1271) 
several times. Setting to zero the coefficients of A, y, A 2 , y 2 and Xy in polynomial S we 
obtain the set of conditions: 

» = -0 + S + iS 
o = --i + S + S 

0 = -1 + 
u 2 ' a -1 

0 = “1+6=7 

0 = ~p + + iri ■ 

The solution is given by the specified values for a, b, c, d and the skew-symmetry condition. 


It was proved in [39] (see also [29]) that under a skew-symmetry condition the obliquely 
reflected Brownian motion does not reach the non-smooth part of the boundary. In the same 
way the above skew-symmetry equality is reminiscent of the second condition in Corollary 
m which is an inequality. We therefore get a partial but handy statement. 


Corollary 23. Assume 


max{a, <5} > ^ 
aS — /Ly > 0 
ay — 7 $ > 0 

lP = Xi- 


and 59 — f3y > 0 


Then the unique solution to [28\) has an invariant distribution given by 


T I 1 + 2a,2a 


59 — (3y 
a5 — j3y 


T 1 + 2(5,25 


ay — 7 6 
a5 — f3y 


Proof. This is a straightforward consequence of Corollary flOl Theorems [TT1 Il3l Il4l and 
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Remark 24. When p = 0, the skew-symmetry condition a/3 + 'yd = 0 means that the 
directions of interaction r x and r y are orthogonal. Then the second condition in Theorem 
m means that ( 9,r]) points into the interior of the quadrant designed by r x and r y . This 
condition ensures recurrence of the process, while the first condition is now a consequence of 
the skew-symmetry condition. 

Remark 25. With the same proof, we may check that under the skew-symmetry condition, 
when 9 = rj = 0, the function q(x, y) = x 2a y 2S is a non-integrable invariant density that does 
not depend on the other parameters. 


A simple change of variables (beta-gamma algebra) provides the following result. 
Corollary 26. With the conditions and notations of Theorem \22l the two-dimensional process 

W t := 

Z t 


(29) 


cX t 


cXt+dYt 

:= cX t + dY t 


has B(a,b) ®T(a + b,l) for invariant distribution on [0,1] x [0, oo) 
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